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Recent theoretical work has derived the correct form of the Ginzburg-Landau differential equa- 
tions, for the superconducting order parameter and vector potential, in the presence of a small 
defect. Here, these equations are applied to the case of a single vortex line pinned on such a defect. 
We develop the coupled set of partial differential equations, and show how to derive analytic solu- 
tions for the order parameter and magnetic field perturbations in the region of space near the defect. 
Certain properties of the unperturbed vortex solution are needed to totally specify our result; these 
are evaluated numerically, and compared with those deduced from Clem's approximate solution. 
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I. INTRODUCTION 

A vortex line in a type-II superconductor induces spa- 
tial variations of the order parameter, ?y(R), and mag- 
netic field, B(R). These quantities can be studied, for 
temperatures close to T c (the superconducting transi- 
tion temperature), and on length scales greater than £ 
(the zero temperature correlation length), by using the 
Ginzburg-Landau (GL) theory. When a single small de- 
fect or impurity is present, the vortex line can lower its 
free energy by positioning its core at the impurity site; 
this effect is known as vorteXj-ninning. nn 

Recent work by Thunebergjm and others,H'Q has shown 
how to correctly include into the GL framework a local- 
ized impurity potential which is not necessaxily weak. 
The pinning energyBtil and pinning potential!!] have al- 
ready been calculated in several situations; such quanti- 
ties are the principle ingredients of macroscopic pinning 
theories. However, the impurity will also influence the 
superconductor in other interesting ways, which can be 
observed experimentally, using microscopic probes. For 
instance, both the order parameter, f?(R), and magnetic 
field, B(R), are altered near the impurity. These per- 
turbations therefore provide important insight into the 
nature of pinning. From a theoretical perspective, the 
calculation of the perturbations in ?7(R) and B(R) con- 
stitutes a non-trivial application of the Thuneberg the- 
ory. 

In this work, we investigate the changes of the order 
parameter 5n(R), and magnetic field <5B(R), due to the 
impurity, when a vortex core is placed on a defect. Near 
the defect these quantities can be calculated analytically. 
Away from the defect, 5r/(R.) decays exponentially over 
the correlation length scale, £(T), while <5B(R) decays 
exponentially over the penetration length scale A(T). 

The paper is organized as follows. In Section II we 
review the appropriate GL theory, and set up the gen- 
eral equations for our problem. These equations, (Elf) 
and (^2|) , are the main results of this section. They are 
a pair of coupled, linear, partial differential equations for 



the impurity-induced changes in the order parameter and 
vector potential. 

In Section III we consider in particular the spatial re- 
gion close to the impurity, where 5r)(R.) and <5B(R) are 
largest. In this region we are able to derive explicit ex- 
pressions for the perturbations — the simplicity of the re- 
sult for 6B(R) is somewhat remarkable. Eqs. ( |25f]27| ) 
contain our main results. In Section IV we discuss the 
results, and estimate the magnitude of the perturbation 
effects. In the appendix we calculate, numerically, the 
GL solutions for vortices in the absence of impurities. 



II. GL EQUATIONS 

A. General Formulation 

For simplicity, in this paper we consider a supercon- 
ductor with a spin singlet, isotropic (i.e. k-independent) 
energy gap. The GL free energy, in terms of the com- 
plex order parameter r/(R) and vector potential A(R), is 
given, in the absence of the impurity, by 
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Here, T> = Vr + 2ieA/hc is the gauge- invariant deriva- 
tive, and the local magnetic field is given by B = Vr x A. 
The coefficients K, a, and j3 can be evaluated by using 
the microscopic theory of a superconducting Fermi liquid. 
The results are 
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where N(0) is the density of states at the Fermi surface, 
and Vp represents the Fermi surface average of the square 
of the Fermi velocity. Note that K can be more generally 
represented as a tensor Kij . However, we have considered 
here the isotropic case, = KSij. 

It is convenient at this point to rescale our variables in 
order to simplify the following equations. We will mea- 
sure lengths in terms of £(T) (with £ 2 (T) = K/2\a\), and 
will normalize the order parameter in terms of its bulk 
value, T)b(T) (with rj 2 = \a\/2j3), for a uniform super- 
conductor. Thus we define the following dimensionless 
quantities: 
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= V r x a = 2eZ 2 (T)B(r)/hc. 
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When an impurity is located at r = 0, the. (dimension- 
less) GL equation acquires an extra term:E§a 

w, _ \ih\ 2 ip+D 2 iP = , aTt \ V X, r. , 

V m W v 576(k B T c ) 3 ^(T)\a\ 

x (D ^ (r)| r =o) • (D <5 3 (r)) . (10) 

To understand this equation, several points should be 
borne in mind: 

(1) The effect of an impurity on the normal state quasi- 
particles can be represented by a potential v(k, k') for 
scattering from k to k' on the Fermi surface. For sim- 
plicity, we take the potential to be of s-wave form, so that 
f(k, k') = v. The impurity parameter er, which, appears 
in the right hand side of Eq. ( |l0| ) , is defined asa 

_ N 2 (0)tt 2 v 2 
° ~ l + 7V 2 (0)7r 2 w 2 ' 

The quasiparticle scattering cross-sectiorusJjhen propor- 
tional to a jk 2 F . In the present theory,!!] oH a "small" 
impurity implies that a/kp -C £q. Thus, the small pa- 
rameter for the expansions which follow becomes a/kp^. 
However there is no restriction on the magnitude of v. 
Note that we consider here a general Fermi surface, sub- 
ject only to the symmetry constraint = KSij ■ For an 
s-wave order parameter and for a Fermi surface with at 
least orthorhombic symmetry, Thuneberg has shown how 
to go beyond the assumption that the impurity potential 
satisfies w(k,k') = vn With Thuneberg's generalization, 
Eq. (|lTJ) remains valid, provided that Eq. ([y]) is gener- 
alized appropriately. 

(2) To derive the GL equation ©, we coarse gram 
a more microscopic theory,El and so lose information on 
length scales shorter than £o = Uvf /ksTc. At the GL 
level, the small impurity appears as a ^-function driving 
term. 



(3) Eq. ( |lC)| ) should be expanded to first order in our 
small parameter cr/kp^ 2 . This amounts to the following 
linearization procedure. We write ip(r) = V>o( r ) + 5ij)(r), 
a(r) = a (r) + <5a(r), and D = V r + ia , where ipo(r) 
and ao(r) solve Eq. ( fL0| ) in the absence of the impurity 
(i.e. with a — 0). We then expand the left hand side 
of (|l0|) to first order in Sip and Sa, and set this first or- 
der piece equal to the right hand side. Sip and 5a. then 
represent the parts of the GL solution proportional to 
a. Consistent with this linearization procedure, the right 
hand side of Eq. (|lfj) involves only the unperturbed terms 
V>o(r) and a (r). 

In the presence of the impurity, the (dimensionless) 
Maxwell equation also picks up an additional term (J 2 ): 

V r x V r x a = 47r(Ji + J 2 ), (12) 

ji= "^S (w_w) ' (i3) 
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Equation ( |l2| ) should also be linearized in Sip(j) and 
<5a(r), with Ji evaluated to first order in these quanti- 
ties and J 2 evaluated using ipo(r) and ao(r). 

For a given (ip (r), a (r)), Eqs. (|ic| ) and (|l2] ) consti- 
tute a pair of linear, coupled, differential equations for 
Sip and (5a. For an application of these equations in a 
simpler context, see Ref. |?]. 

B. Application to a Single Vortex Line 

We now specialize, and apply Eqs. ( |l0| ) and ( |i~2| ) to 
the case of a single vortex line which is centered on an 
impurity at the origin. ipo(r) and ao(r) then represent 
the solutions of the a = GL equations, with a single 
vortex at the origin. These unperturbed vortex solutions 
are well-studied.El (Some aspects of the computations are 
discussed in the appendix.) The following gauge choice 
is most convenient: 

Mr) = Me-**, (15) 
ao(r) = a (p)<i>, (16) 

where we have used polar coordinates (a, (p, z) and taken 
the vortex line to be along the z-axis.Q At large p, we 
have fo(p) — > 1, while in the limit p — > 0, we have 

lim{/o(p) = 7P, a (p) = rp). (17) 

7 and r are constants which depend on (3 and K through 
the Ginzburg-Landau parameter n 2 — A 2 (T)/^ 2 (T) = 
Ti 2 c 2 f3/Aire 2 K 2 . (7 and r should be determined numeri- 
cally, as described in the appendix.) 

According to Eqs. ( |To| ) and (JTJ), the impurity pertur- 
bations to ipo (r) and ao(r) take on the following form: 



2 



61>{T) = f 1 {p,z)e- i *, 
<5a(r) = ai(p,z)(f>. 



(18) 
(19) 



Thus, the order parameter magnitude will depend on z 
as well as p. Further, we find 



<5b(r) = V r x 
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-— p + --(, ai) z. (20) 



So, the magnetic field in the presence of the impurity will 
now have a p component, in addition to a perturbed z 
component. 

Using the preceding definitions (|l5|), (|l6|), ([l8|), and 
(|l9|), we arrive at the following pair of coupled, linear, 
differential equations for fi(p, z) and a±(p, z): 
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Note that the magnitudes of /i and eti are determined 
by the (^-function driving term in Eq. (pl|), particularly 
through the impurity parameter a. We also note that 
since the unperturbed vortex line has a vanishing current 
density at the impurity site, the explicit impurity term 
( |l4| ) in Maxwell's equation is zero. Thus, the changes in 
the magnetic field are due solely to the changes in the 
order parameter caused by the impurity. 



III. SOLUTION CLOSE TO THE IMPURITY 

The rather formidable set of equations for fi(p, z) and 
ai(p, z) simplifies greatly, at distances much closer to the 
impurity than £(T); in our rescaled variables, this means 
that the results of this section are meant to cover the 
region /> 2 + z 2 «l. Recalling that fo(p) and clq(p) both 
vanish linearly as p — > 0, the proper short distance ver- 
sions of Eqs. (|l]) and ( |22| ) become 



d 2 h , 1 dh , 8 2 h h 



+ 



dp 2 p dp dz 2 p 2 



576(M) 3 £ 5 (T)H 



r (x-iy)- V r 5 3 (r), (23) 



d 2 ai 1 da\ d 2 ai 

H H K~ 

dp 2 p dp dz 2 



— = -—fx 

2 tc 1 



p 



(24) 



One advantage of the formulation we have chosen is 
that the equations become decoupled in the following 



sense: we can now solve Eq. ( J23| ) independently for fi, 
and then use this result in Eq. (pi]) to find ai. Re- 
markably enough, when the analytic form for fi(p,z) 
is inserted into Eq. (|24|), a simple, analytic solution for 
ai(p,z) may also be found, leading to a simple solution 
for <5b(r). The following results are obtained: 
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IV. DISCUSSION 

Eqs. (p5|-p7[) are the main results of this paper. They 
reflect the changes in the order parameter, magnetic field, 
and vector potential due to the impurity. We now discuss 
important points concerning these equations: 

(1) The results hold for distances from the impurity 
greater than £ , but less than £(T). Close enough to 
T c , this range of validity can be reasonably large, since 
£(T) - f [(T c -T)/T c ]- x / 2 . Furthermore, any anomalies 
which appear in Eqs. ( f25f|27| ) as p, z — > are not physical, 
since these solutions are not valid at short distances. 

(2) The magnitude of the order parameter near the vor- 
tex core actually increases from its unperturbed value, 
when the impurity is taken into account. This is con- 
sistent with the idea that scattering from nonmagnetic 
impurities lowers_the free energy cost of gradients of the 
order parameter a The situation is analagous to the case 
of a finite concentration of impurities, where the GL co- 
efficient K is reduced from its impurity-free value. 

(3) Associated with the increase in |^|, circulating vor- 
tex currents are also enhanced near the defect site, caus- 
ing the z component of b(r) to increase. We may esti- 
mate the largest deviation of the magnetic field from its 
unperturbed value, <5-B max - To do this, we estimate <5b at 
a distance £o from the impurity. Note that in real units, 
B = H c2 {T)b, where H c2 (T) = hc/2e£ 2 {T). If we take 
7 2 and a to be of order one, we finally get 
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It should be possible to test both the magnitude and tem- 
perature dependence of this prediction using a suitable 
microscopic probe£3 

(4) In spite of the local perturbations to the magnetic 
field, impurities do not affect the quantization of vortic- 
ity; the net flux associated with a pinned vortex remains 
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<f>o in the z direction. However, the vortex field lines, 
which are parallel to z for the unperturbed vortex, now 
acquire a new p component near the defect. The pertur- 
bation to the magnetic field is largest inside a radius £o 
of the defect [see Eq. (p8|)] . At a distance of order £ (T) 
from the vortex core, the perturbation field lines begin 
to close on themselves. Full screening of <$B(r) occurs at 
distances of order A(T). Note that neither the screening 
of £B(r), nor the decay of 5tp(r) over length scale £(T), 
can emerge from Eqs. ( |23| ) and (24), since those apply 
only very near the impurity. 
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APPENDIX: DETERMINATION OF 7 AND r 
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FIG. 1. Results are shown for 7 (bottom two curves) and 
t (top three curves), denned in Eq. (0), as a function of the 
Ginzburg-Landau parameter k. Exact numerical solutions of 
Eqs. (Al) and (A2) are given as solid lines. Approximate solu- 
tions, due to Clem, are shown as dashed lines; these describe 
t very accurately for all n. The large-K asymptotic form of 
Clem's result, r ~ (\xi\/2k — 7_b)/2k 2 , is shown as a dotted 
line, and is accurate in the range k •> 5. 



In Section IV of this paper, an estimate was made for 
the magnitude of 5B, near the defect site, for which it 
was found that 5B oc 7 2 . In general, many quantities 
calculated near the vortex core depend sensitively on the 
parameter 7 = dfo/dp\ p ^,o, regardless of the presence of 
a perturbing defect. To the best of our knowledge, there 
exists in the literature no general calculation of this im- 
portant quantity, nor of the quantity r = da®/ '<9p|p_>o- 
Such a calculation is therefore presented here. The re- 
sults are found to_agree quite well with the approximate 
solution of Clem.Lj 

Using dimensionless variables, the equations satisfied 
by fo(p) and ao(p) become 
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while the boundary conditions, appropriate for a vortex 
core located at p = 0, are fo(fl) — ao(p)p — when 
p — ► 0, and fo(p) — ao{p)p — 1 when p — » 00. Eqs. (Al) 
and (A2) can not be solved exactly, in the general case. 
We therefore solve the equations numerically, for fixed 
values of k. To characterize our results, we show the re- 
sulting k dependences of 7 and r in Fig. 1, with k > 1/V2 
for a type-II superconductor. 



The results of Fig. 1 reflect an exact treatment of 
the vortex, within the GL description. An approximate, 
but analytic, treatment has also be provided by Clem.til 
In this scheme, only the second GL equation (|A2|) is 
solved. This is accomplished by using the variational 
form /o ~ p/(p 2 +(,v) 1 ^ 2 ) which permits an exact solution 
for ao- The GL free energy is then constructed for the 
vortex line, and minimized, in terms of the variational pa- 
rameter £„. Clem's model thus improves and regularizes 
the London vortex solution, by introducing a vortex core. 
Indeed, the results for the quantity r (which emerges 
from the conventional London theory only through ad hoc 
regularization) are quite impressive; Fig. 1 shows Clem's 
approximate, transcendental solution for r, as well as its 
large-K asymptotic behavior: t ~ (lni/2K — 7£j)/2/t 2 , 
where 7^ ~ 0.577 is Euler's constant. Although the 
model is not designed to provide accurate information 
in the vicinity of £(T), we may still use it to calculate 7 
as shown in Fig. 1 . It is interesting that the outcome has 
the correct qualitative behavior, but is slightly too large. 
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